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Abstract: 

Polymer chains with hard-core interaction on a two-dimensional lattice are modeled by 
directed random walks. Two models, one with intersecting walks (IW) and another with 
non-intersecting walks (NIW) are presented, solved and compared. The exact solution of 
the two models, based on a representation using Grassmann variables, leads, surprisingly, to 
the same analytic expression for the polymer density and identical phase diagrams. There 
are three different phases as a function of hopping probability and single site monomer 
occupancy, with a transition from the dense polymer system to a polymer liquid (A) and 
a transition from the liquid to an empty lattice (B). Within the liquid phase there exists a 
self-dual line with peculiar properties. The derivative of polymer density with respect to 
the single site monomer occupancy diverges at transitions A and B, but is smooth across 
and along the self-dual line. The density-density correlation function along the direction 
x, perpendicular to the axis of directedness has a power law decay 1/x 2 in the entire 
liquid phase, in both models. The difference between the two models shows up only in the 
behavior of the correlation function along the self-dual line: it decays exponentially in the 
IW model and as 1/x 4 in the NIW model. 
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A directed random walker makes steps along a given (z) axis only in one (forward) direction. 
Random fluctuations are present in the transverse directions 1 . Such walks are the subject 
of great interest since they may model diverse phenomena like polymers under flow 2 , tracer 
diffusion 3 , electrorheological fluids 4 , commensurate - incommensurate phase transitions 5 , 
vortex lines in high T c superconductors 6 ' 7 , world lines of quantum particles 8 , the behavior 
of interfaces in 1 + 1 dimensions 9 , some aspects of biomembranes 10 , etc. They share a 
number of properties of dimers 9 and vortex models 11 . 

Directedness makes analytic calculations considerably easier, and a number of exact 
results exist both for a single walk and cases with many walks. 

In the present work we concentrate on some generic properties of random walks. 
We introduce and solve exactly a model with non-intersecting walks. Using the same 
formalism, we also obtain the solution of a similar model, which allows for intersections. 
The comparison of the two models reveals some surprising effects. The interesting feature 
of our systems is that along with second order phase transitions, observed in previous 
investigations of similar models 12-14 , they contain a self-dual line, with characteristics of 
yet another critical phenomenon. 

We first discuss the NIW model depicted in Fig.l. Chains are pinned with one of their 
ends at z = 0. In order to avoid crossing of walks in the vertical direction only steps of 
length 2a are allowed, where a is the lattice spacing. Each such chain can be parametrized 
with two indices, r and j. Here r are lattice sites denoted by open circles in Fig.l. The 
index j takes values or 1. j = if a site is an r-site, and, as a consequence, can be 
reached from another r-site making two steps along the elementary vectors e\ (site 0\ in 
Fig.l) or e 3 (site 2 in Fig.l). j = 1 for a site which can be reached from an r-site by 
making a single step in the direction of the vector es (site O3 in Fig.l) or from a j = 1 
site directly below (site O4 in Fig.l). Chains start at z = 0, at r-sites (also needed to 
avoid crossing). This is model A of ref. 14 which was used to describe phase transitions 
in biomembranes. In the present article we investigate correlations in this model (as well 
as in the IW model) to provide a deeper understanding of the statistics than what can be 
obtained from the thermodynamic properties studied in ref. 14. 

We impose periodic boundary conditions both in the x and z directions. We assume 
a hard-core potential between walks: chains repel each other with infinite energy upon 
contact. Therefore, the contribution of the configuration in Fig.l to the partition function 
is zero because of point O4. We assign unit weight to a step in the vertical direction 
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(with length 2a along the vector ei), and a weight w to a step along the diagonals of the 
elementary square. The weight (fugacity) of an empty site is denoted by fx. Introducing 
a pair of Grassmann variables 15 {^r,j, ^r,j} for each lattice site, the configurations of our 
system can be generated from the partition function 

Z = [ eX P (^ ^'rJ'r/r'^'r'.rr^^ = ^tG^ . (1) 

J r,r' j=0,l 

Here T>^>T>^> denotes the Berezin integration 15 over all Grassmannians, and, the lattice 
Green's function is given by its Fourier components 

x-i _ ( e ikl - /j we ikl (l + e ik2 )\ , . 

Gfc - { w (l + e - ik *) e^-iJ J • { ] 

The 2 x 2-matrix structure reflects the dependence on the index j = 0, 1; ki and k 2 are 
the components of the two-dimensional momentum vector along the x- and ^-directions, 
respectively. On a lattice with iV sites expression (2) leads to a free energy 

llogZ = f f f^ff log(^ 2 - 4«,V fc W ^ - 2^ + e 2 ^). (3) 

In what follows we first outline the calculation of the density n r of lines at site r and 
of the density-density correlation function. The probability that a line goes through the 
site (r, j) is given by 

n r = (1 - ^o^V.o) = n, (4) 

where the last equation follows from translational invariance. The average value in (4) is 
to be calculated using the partition function given by (1) and (3). It is easy to see that n 
can be expressed in terms of the Green's function as 

n = 1 + ^G r , r; oo- (5) 

The calculation of the Green's function is tedious but straightforward. Starting from the 
Fourier components of the Green's function in (2), the diagonal matrix elements of G read 

' r;0 ° J- n J- n (n - e^) 2 - 2w 2 e^ (1 + cos k 2 ) 2n 2n ' 1 J 

By evaluating the double integral in (6) we finally arrive at 
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n = 1 ^— cos 1 

2tt 



1 



(1 + ^ 



1 



(1-/*) 
2w 2 



2-, 



(7) 



2(fi + w 2 ) 

Here =F corresponds, respectively, to cases \i > 1 and < 1. 

The result given in (7) corresponds to the phase diagram shown in Fig. 2. \i = 1 ± 2w 
(A and B respectively in Fig. 2) are lines of phase transitions from a dense system of 
directed polymers (n = 1) to a polymer liquid (0 < n < 1) and finally to a system without 
directed polymers. Approaching the phase transition lines from the regions II and III along 
directions parallel to the w-axis, one finds dn/dji ~ l/v^e, where e = 1 — (1 — fi)/2w and 
e = 1 — (fj, — l)/2w, respectively for the two cases. There is another critical phenomenon 
at ji = 1 (broken line in Fig. 2). Across this line the above partial derivative of the 
density varies smoothly. This critical behavior can be understood from the density-density 
correlation function C X2 , which measure the correlation perpendicular to the polymer 
direction. C X2 is given in terms of the Green's function 



(/j - e ikl )e ik2X 



dki dk2 



n (fi-e ik i) 2 -2w 2 e ik i(l + cosk 2 ) 2n 2n 



(8) 



as 



C X = ^G X G_ X . (9) 
Performing the k\ integration in (8), one arrives at 

g x = i - r Mi-»)+* Pik2X dk2 (10) 

2 J-k* A/Q^+4^a 2n 

with k* = cos" 1 } ^ 1 ^ 1] and a = 2w 2 (l + cosk 2 ). k* takes the values 7r and along 

the line = 1 and the transition lines \i = 1 ± 2w, respectively. The calculation of the 
integral in (10) is an exercise in complex contour integration, and in the large x limit can 
be performed analytically for any value of k* , with the result 



(-iy 



for k* = n 



~ < ,* ,u. s ■ (11) 

Finally, the density-density correlation function decays as l/x 4 along the line fx = 1 and 
as 1/x 2 elsewhere in the liquid phase. 

The special property of the line ji = 1 is also elucidated by noting that the expressions 
for the density of polymer lines, given in (7) are invariant under a "duality transformation" 



4 



[i — > 1/n and w ^ w/n. This duality at = 1 reflects the equivalence of directed polymer 
lines and lines of empty sites on the lattice, /j, = 1 is a self-dual line along which polymers 
and empty sites compete with the same weight. This competition leads to a faster decay 
of the density-density correlation function. 

In regions I and IV in Fig. 2, n has the values 1 and 0, respectively. Fig. 3 shows the 
behavior of n along the lines w — 0.4 and w = 0.6. 

If the restrictions imposed on the model as defined above (only steps of length 2a are 
allowed in the vertical direction, chains must start at z = 0, at r-sites) are relaxed, the 
model describes a system of intersecting random walks. The analogous calculations for 
the partition function (Z), density of lines (ft) and the correlation function (C x ) are even 
simpler in this case (one deals with a simple Bravais lattice as opposed to a lattice with 
basis in the case of non-intersecting walks). We only quote the results 



Although the above expression for the free energy clearly differs from (1), surprisingly, for 
the density one obtains ft = n, where the expression for n is given in (7). This in turn 
implies that the phase diagram for the IW model is the sane as the one obtained for the 
NIW model. Finally, the correlation function C x can be calculated similarly to C x in (9). 
However, instead of G x in (10) we now have for \i = 1 



with a = 2wcos/c2- Considering /i = 1 and, for simplicity, w = 1 and comparing G, 



different. As a function of z = e lk2 the square root of G x is analytic in a vicinity of the 
unit circle. Therefore, the path of integration can be contracted to a smaller circle with 
\z\ < 1. Consequently, G x decays exponentially on the \x = 1 line. Everywhere else in the 
liquid phase 13 it decays like 1/x. 

In conclusion, the analysis of the simple models of non-intersecting and intersecting 
directed polymer chains presented in this work (and studied earlier in references 13 and 14) 
leads, surprisingly, to identical thermodynamic properties. The phase diagram is unusually 
rich with multiple transition lines and phases. There are second order phase transitions, 
with diverging second derivatives of the free energy. The region of the phase diagram 
between these transitions, the liquid phase (A and B in Fig. 2) is critical in the sense that 




(12) 




(13) 



and G x we notice that the arguments of the square root in the denominators are quite 
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in it the density- density correlation function decays everywhere according to a power law. 
An additional feature of the phase diagram is a self-dual line in the liquid phase, which 
manifests itself also in the change of the decay of correlations. This change of correlations 
does not seem to affect the behavior of the thermodynamic quantities; they vary smoothly 
across the self-dual line. The difference between the two models, in our study shows up 
only in the behavior of the correlation function along this line. 
One of us (G.F.) acknowledges useful discussions with S. Bhattacharjee. 
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Figure Captions 



Fig. 1 The schematic representation of the lattice where the model presented in this work 
is defined. The figure shows a typical configuration of some of the directed polymers. The 
use of periodic boundary conditions in the z direction is indicated explicitly. For more 
details see text. 

Fig. 2. The phase diagram of the model. In regions I and IV, n = 1 and n = 0, respectively. 

Fig. 3. The behavior of the density of polymer lines along lines w = 0.4 (full line) and 
w = 0.6 (dashed line) in Fig. 2. 
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